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ABSIRACT 

This report describes the research conducted over the past seven years by three groups at 
the University of California, San Diego, on the seismotecionics of plate boundaries. The pro- 
ject involved both instrumental development and an observational program designed to study 
various aspects of the seismotectonies of southern California and the northern Gulf of Califor- 
nia Dr Jon Berger and Frank Wyatt were primarily involved with the development of the 
Cecil and Ida Green Pifion Flat Observatory. Dr James N. Brune and his group worked on the 
interpretation of the seismotectonies of southern California and Baja California. Dr. John 
(lotxlkind and Dr. Richard Warburton further developed a unique su{>ereonducting gravimeter 
and were supported under this program for the deployment and operation of these instruments 
at several sites. Each of these three sub-programs are reviewed separately in Sections 2, 3. and 
4 of this report while Section > discusses the work of Dr. Duncan .Agnew and Dr. Christopher 
Beaumont on earth tides. 
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I.O Introduction to the Report 

This project, involving three groups at the University of California, San Diego, was basi- 
cally an instrumental and observational program designed to increase our knowledge of the 
seismotectonics of plate boundaries. Dr. Jon Berger and Mr. Frank Wyatt were primarily 
involved with the development of the Cecil and Ida Green Pifton Flat Observatory (PTO). Dr. 
James N. Brune and his group worked on the interpretation of the seismotectonics of southern 
California and Bjya California. Dr. John Goodkind and Dr, Richard Warburton further 
develoned a unique superconducting gravimeter and were supported under this program for the 
deployment and operation of these instruments at several sites. Each of these three sub- 
programs will be reviewed separately in Sections 2, 3, and 4 of this report while Section 5 
discusses the work of Dr. Duncan Agnew and Dr. Chris Beaumont on earth tides. 

Each of these projects was, in different ways, involved in the support of NASA’s study of 
the seismotectonics of the San Andreas Fault System. 

The general teC''n' :s of the San Andreas Fault (SAF) system, and the basic mechanics of 
faulting and earthq .s, have been fairly well established in recent years. The SAF forms part 
of the boundary beiween the Pacific plate and the North American plate. Faulting o;curs as a 
result of shearing between these plates as they drift past one another upon a weak, lubricating 
asthenosphere at a depth of about 50-100 km The frictional resistance resulting from the pres- 
sure of one side of the fault against the other causes the fault to "chatter" on a large scale. 
Str.iin energy is stored while the fault is locked and is icleascd when slip occurs. 

Slip along the fault which relieves these stresses can occur seismicaily or aseismically. 
Slow slip (creep) does not produce earthquakes under normal circumstances. However, if the 
slip prtxeeds at a speed comparable with the velocity of propagation of elastic waves then 
seismic waves are excited and an earthquake occurs. Since earthquakes are highly dangerous to 
life and property, it is essential to understand the underlying differences between these two 
mtxles of strain release. The estimation of seismic risk, the prediction of earthquakes and ulti- 
mately the prevention of disastrous events all hinge on a thorough understanding of fault 
motion mechanisms. Specifically, we must determine the processes by which: 1) strain energy 
IS stored in the earth preceding large earthquakes, and 2) strain energy is released by creep. 

Our studies have involved monitoring strain accumulation of some southern California 
faults which form at the plate boundaries. Here we differentiate between 1) the high seismic 
risk surficial creep, or creep at depth where the surficial region is locked; and 2) the low seismic 
risk creep that occurs at low stress throughout the plate boundary. 

The major uncertainties that remain in our mechanical picture of the plate boundary are: 

On a large scale. 

1 What arc the driving forces responsible for the plate motions and how are they dis- 
tributed in spiice and time'!' 

2 . What are the current rates of relative plate motion both across San Andreas Fault 
(SAF) and at the centers of the plates? Do the motions agree with one another and 
with the geological rate of 5-6 cm/yi 

3 What IS the tectonic relationship between fault motion along the SAF in California 
and the opening of the Gulf of California'.’ Is this motion occurring predominantly 
seismicaily or aseismically'’ 

.And on a smaller scale. 

4 If there is a significant difference between the motions of the plate centers and the 
motions at the boundary, where docs the plate distortion occur? Is it near the fault, 
distributed throughout the plate, or on the fault rones'.’ 

5 Mow IS aseismic slip accomplished below the /one of earthquakes’ Is it steady or 
episodic ’ Does a significant part of the energy in large earthquakes come from this 
/one'’ Are large earthquakes triggered by premonitory deformation in this /one ’ 
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6 Is there a significant amount of aseismic movement aa'ompanying earthquakes and 
if so how can we measure the total moment? 

7 Are there premonitory changes in the elastic properties in earthquak source regions 
as the dilatancy models suggest, and can these changes be used to predict earth- 
quakes? 

Recent development of theories of plate tectonics point to the Gulf of California as one cf 
the most important areas in which to study the mechanism and characteristics of plate boun- 
daries and ocean-floor spreading process The Gulf is underlain by an oceanic ridge and is itself 
the result of sea-floor spreading processes related to this ridge. Nowhere else in the world is 
such an area so easily accessible to study with detailed land-based seismic and geodetic arrays. 
At the same time, the Gulfs evolution has progressed to the point where its spreading 
prtK-esses are probably very similar to those of ridges in the centers of the rntyor ocean basins. 
The tectonic history of the Gulf of California is intimately related to that of the San Andreas 
fault system and full understanding of this system demands a better knowledge of the mechan- 
ics of the ocean-floor processes at its two extremities. 

In the sections that follow, we will describe the activities supported by the grant under the 
headings; 

Section 2. The seismotectonics of southern California — Pifton Flat Observatory. 

Section 3. Seismicity studies related to the Gulf of California. 

Section 4. Development of the superconducting gravimeter. 

Section 5. Earth tide studies. 
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2.0 The Seismotectonics of Southern CallfomU: PiAon Flat Observatory 

There are several distinct methods of investigating the seismo-tectonics of the earth’s 
crust. Traditionally, seismic studies have perhaps yielded the most results. Here, naturally 
occurring earthquakes are the sources used to investigate both the earthquake process itself and 
the genera! tectonic fabric. Recently, geodetic surveys have been used rather extensively in 
California to study the spatial and temporal pattern of strain accumulation in the crust, If the 
crust is an elastic material, the strain field will be simply related to the stress field and that in 
turn is a key factor in understanding the plate motions in general and the earthquake process in 
particular. A third method of investigation is the development of continuous rec'>rding instru- 
mentation of high sensitivity and stability to measure the surface deformations directly. 

This element of the NASA supported program involved the development and operation of 
precision instrumentation to monitor crustal deformations. For Dr. Berger’s group, it 
represented a continuation of work begun at Camp Elliott, near La Jolla, in the late sixties on 
the laser strainmeter (Berger and Lovberg, 1969). In 1971, we began *o deploy this type of 
instrumentation in a more geophysically interesting site with three m^or goals in mind: 

1 . The development of new instrumentation to measure crustal dynamics. 

2. The establishment of a reterence station consisting of the most stable and sensitive 
instrumentation available in an attempt to determine true ground deformation rates 
and thus improve the limits of detectability of anomalous tectonic behavior. 

3. The monitoring of the tectonic and seismic activity near PFO — an area of unusually 
high seismic activity. 

These goals are quite fundamental to our understanding of th? dynamics of the faulting 
process. Further, the specific program embarked upon and supported by NASA is a necessary 
element in the interpretation of geodetic measurements of crustal deformation. While the scale 
of the measurements (—1 km) is small by geodetic standards, the problems addressed are com- 
mon to any observations of crustal dynamics. The goal of monitoring stress changes in the 
earth's crust is currently limited by the noise in the observations--noise which tends to mask 
tectonic and seismic effects. Our efforts are directed towards identifying, reducing and, when 
possible, eliminating some of these noise sources in order to examine the underlying signals of 
interest. 

It is perhaps worth reviewing the relevance of this research to NASA’s goals in Crustal 
Dynamics and related studies. NASA's involvement in this subject arises from the techniques 
they have developed to measure the position of points on the earth’s surface relative tc earth 
satellites and extragalactic sources. The research goals of these studies relate to seismic and 
tectonic activity, polar motion and rotation rate. 

Our goals at PFO are directed towards the understanding of both seismic and tectonic 
motions and their interconnections. We are pursuing an observation program that involves the 
development and testing of new instrumentation and measurement techniques. We are observ- 
ing strain, tilt and gravity continuously with what we believe are the most sensitive and stable 
instruments currently available. This program has, over the past few years, provided consistent 
upper limits to the spectrum of crustal deformations occurring (at least) in the area of the 
observatory. This is an area of high tectonic activity, high seismic risk, and is very likely typical 
of southern California 

While our measurements are made over rather short (by NASA standards) baselines, the 
sensitivity and stability is rather high Figure 2.1 compares the capabilities of observatory meas- 
urements, land based EDM (Electromagnetic Distance Measurement) techniques (such as util- 
ized by the USGS m their trilateration networks), and extraterrestrial surveying techniques. 
The trade off between baselength and resolution is obvious. But what is not obvious is the spa- 
tial aliasing and temporal aliasing to which long baseline, infrequent measurements might be 
subject. 

In the case of the USGS trilateration networks (see Savage ct al.. 1981). the survey areas 
are irregular figures with typical dimensions of 100 km (The Palmdale net is a notable 
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Figure 2 \ A comparison of siram measuring techniques 
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exception.) The time taken to survey a single network comptitely is usually many weeks. The 
data are reduced under the assumption of spatially uniform strain and constant strain rate 
(Prescott cf ai, 1979). Whether the tectonic strain fields approximate the assumptions of spa- 
tial uniformity and temporal rate constancy has yet to be detcmined, but the impact of these 
assumptions to NASA’s Crustal Dynamic; observational program is considerable. 


What strain rates should we expect to see in such a tectonically active area as southern 
California? The rather substantial body of trilateration data collected in this area over the past 
seven years or so (Savage et ai, 1981) clearly indicates that the normal secular strain rates 
(normal over this period at least) are 1-2 x lO Vyr. The NASA SAFE experiment monitoring 
the length of a line from Otay Mountain near San Diego to Quincy, California observed on the 
average a strain rate of lO’^/yr over that 900 km distance (Smith, 1980). Our observations 
indicate that on a 1 km baiielength the same result is obtained. Crudely, then, the strain rates 
of about 10 Vyr are observed over baselines of 1 km to 1000 km at least in this area. 


Ructuations about these average strain rates might be expected to occur if the tectonic 
stress is not accumulating uniformly in time. Obviously, the occurrence of earthquakes intro- 
duce sudden perturbations into the ambient stress field if only due to the local release of strain 
energy. If an earthquake occurs, roughly speaking, the strain change will be proportional to the 
stress drop. 


A( 
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is a constant close to one 
is the stress drop, typically, 10 bars (10^ N/m^) 
is the rigidity usually taken as 3 x 10'° N/m^ 
is the earthquake rupture dimension 
is the hypocentral distance 


As the stress drop is independent of seismic moment, bigger earthquakes do not neces- 
sarily involve bigger strain changes; only the area affected is bigger. As an example, for a 
moderate earthquake in California such as the Borrego Mountain, recent Imperial Valley, or 
San Fernando earthquakes, th** pertinent dimension is about 10 km to 15 km, the depth of 
faulting. Farther away than this, the strains will be smaller than 10"^. Figure 2.2 shows the 
expected average strains as a function of distance and seismic moment (and magnitude) 


We might expect that coseismic strain changes and post seismic strain changes would be 
on the same order of magnitude as a result of redistribution of stress in the lower lithosphere. 
Post seismic surface slip about the same size at the co.seismic slip has been reported (Crook et 
al. 1980) for the 1979 Imperial Valley earthquake, and strain observations strongly suggest an 
exponential strain change with a time constant of about 14 hrs and an amplitude of ~.3 times 
the coseismic strain change for that event (Wyatt, 1982b). Preseismic strain changes, if they 
occur, will not likely be much larger. Thus, we conclude that unless instruments are located 
within the near field of an earthquake, strains of more than a few times 10"^ should not be 
expected. Our observational and instrumental strategies should be based on this premise. 


The results of this premise are: 


1 Extraterrestrial geodetic techniques with an accuracy goal of 1 cm will likely produce 
useful measurements of the strain field related to the faulting process only over 
baselines of 100 km and greater. This applies to existing Very Long Baseline Inter- 
ferometry (VLBl) and Satellite Laser Ranging (SLR) techniques as well as the pro- 
posed Global Positioning Systems (GPS) 

2. Terrestrial EDM (Electronic Distance Measurement) techniques currently limited to 
statistical errors of a few parts in 10 ’ (in the case of trilateration nets or two color 
ranging devices) must be improved in order to resolve strains at least an order nf 
magnitude smaller. Further, it is possible that the spatial averaging and infrequent 
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sampling of the irilateratiori nets may obscure many tectonic signals. If this is so, 
then densifkation and more frequent re-surveys of the networks are needed. 

3 Observatory based instrumentation must be improved to provide lower noise leveis- 
•primarily to distinguish environmental signals from tectonic signals. 

2.1 Seismotectonic Setting 

Pinyun Rat lies at the junction of the San Jacinto Mountains and the Santa Rosa Moun- 
tains, in the northern section of the penin. ular ranges of southern California. Located some 20 
km south of Palm Springs, California (Figure 2.3), it is between the San Andieas and the San 
Jacinto fault zones in an area of considerable tectonic interest. 

The behavior of the aseismic stretch of the San Andreas Fault south of F^rkfleld is of 
considerable importance in view of its past seismic history (the 1857 mag > 8 Fort Tejon earth- 
quake) and the high density of population in the area of southern California. It has long been 
recognized that somehow strains must oe cccumulating in this general area to accommodate the 
inter-plate motions that are so clearly evident both to the northwest and southea.st. The trend 
of the San Andreas Fault swings from NW-SE to nearly EW as it cuts through the southern end 
of the Transverse Mountain Range. This change of direction may present an impediment to 
the gradual release of strain accumulation by creep and small seismic events. 

The San Jacinto fault, however, offers a contemporary alternate fault which can accommo- 
date at least some of the strain accumulation. Indeed, in many respects, this fault seems a 
more likely candidate for the principal plate boundary than the San Andreas. As Sharp (1967) 
has pointed out. the San Jaemto fault is the only fault that cuts the Transverse Ranges and yet 
maintains the northwest-southwest trend that the plate boundary has on the large scale. 
Further, it is clearly the most seismically active fault in southern California, having been rup- 
tured by eight moderate (6 < mag < 7) earthquakes since 1897. (Allen ?f al., 1965; Thatcher 
et al., 19''5>. In contrast, the San Andreas fault has experienced only one mag > 6 event 
south of the great 1857 earthquake in historic times: the 1948 Desen Ho'. Springs event. 

The problem, however, with assigning the San Jacinto Fault as the ntqjor plate boundary 
is the magnitude of the cumulative horizontal motion which has taken place across it. Sharp 
(1967 and 1981) has estimated the total slip in this fault at 24 km of right lateral motion. Pleis- 
tocene slip rates are estimated at 0,8 to 1.2 cm/yr, considerably less than the 6 cni/yr deduced 
as the average North American-Pacific Plate motion in this area (Minster and Jordan. 1978). 
However, Sharp (1980) also points out that in the last few hundred years, there is evidence 
that the slip rate on the San Jacinto Fault may have been as great as 5 cm/yr. He concludes 
that. "If the San Jacinto and adjacent segments of the San Andreas fault alternately assume 
dominant roles in absorbing motion between the Pacific and American plates, perhaps even 
more recently than 400 years ago. the San Andreas fault in the northern Salton Trough 
expres.scd most of the motion but has since become relatively inactive." 

In their study of the historic sei^.micity along the San Jacinto Fault, Thatcher et al. (1975) 
demonstrated that sequence of events along the fault since 1890 has ruptured the fault along 
nearly its entire length. Two gaps apparently remain, one at the northern end of the fault and 
the other in the region of the town of Anza (Figure 1.3) There was a moderate earthquake 
near Anza in 1937. but its moment of 2 x 10’^ dync-cm. indicates that it could hardly be con 
sidered a major rupture of the area. The U S. Geological Survey periodically rc-surveys some 
seven trilateration networks in southern California utilizing laser distance measuring devices 
The networks are irregular figures with typical dimensions of 100 km and it usually takes many 
weeks to complete one resurvey. For each network, the individual line length changes arc least 
squares fit to a model of homogeneous strain where it is assumed that the strain field is uni- 
form over the network and over the interval (i e . in space and time) (Prescott et al., 1979) 

Savage .'/ al., 1981, reported their most recent results which span the years 1973-1980 
Figure 2 4 shows their results for the principal strain rales averaged over this period. The aver- 
age strain rates for all of southern California dunng this period were 18 x 10 Vyr. north-south 
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Figure 2 4 Souihem Catifomia geodimcter nets 
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'.'ontraciion coupled with a 0.8 x lO'Vyr. east-west extension. Howeve'', Savage e( ai, 1981, 
point out that certain regional variations in this picture are suggested by the data. Prior to the 
1978-79 interval, the strain rates averaged to 2.6 x 10~* north-south contraction with an 
undetectable east-west extension. The change in the strain rates appeared first on the southern 
networks and then migrated northwards. It is interesting to note that the seismicity of the area 
for several years prior to the 1978-1979 period was anomalously low (Bufe and Toppozada, 
1980) and since then has been significantly higher. 

With the tectonic and seismic setting in mind, we chose the exact site of the observatory 
(PFO) from more mundane considerations. We sought a suitable site for the installation of 
long baseline (1 km) instrumentation which translated into a requirement for a fairly flat piece 
of land of approximately 160 acres. Further, it was important to locate the observatory on com- 
petent material away from large accumi iations of alluvium. Finally, we sought a relatively 
remote area to reduce cultural interference, and to make the acquisition of suitably sized parcel 
of land a practical matter. 

Pinyon Rat, which is located some 12 km northeast of the San Jacinto fault system and 
25 km southwest of the San Andreas fauk zone (33.609°N, 116.455*W), is an area that has the 
desired properties. The flat is situated on the northern-most flank of the Santa Rosa Mountains 
in the Peninsular Range Province. Roughly 12 km^ in area, the flat is bounded by Asbestos 
Mountain to the north. Deep Canyon to the east, Santa Rosa Mountain to the south, and Palm 
Canyon to the west. Figure 2.5 (after Wyatt, 1982a) shows a generalized map of the pertinent 
geological features and a sectional profile of the area, indicating the smooth topography of the 
flat between "PFO" and "Edge." 

As Wyatt (1982a) points out, there are several possible explanations as to why such an 
isolated flat surface is found at these elevations (1300 m). Undoubtedly part of the reason is 
the nature of the crystalline rocks which underlie the area. The rocks are composed of Pinyon 
Rat Granodiorite (Foster, 1976) of middle Cretaceous age (90-100 Ma b.p.). This rock is 
unusually rich in biotite which causes rapid weathering because of its expansion when exposed 
to moisture (liherwood, 1976). Weathering causes the near surface rock to decompose (via 
granular disintegration) to grus. 

Jahns (1954) suggested three possible mechanisms for the development of low relief sur- 
faces on pre-Cenozoic terrain of crystalline rocks in the Peninsular Range Province: 1) The 
areas may be remnants of a single extensive peneplain, called the Southern California 
peneplain, that have been differentially uplifted by Quaternary faulting. 2) Some of the sur- 
faces may be exhumed features that were originally formed in pre-Quaternary time and only 
recently re-exposed. Some support for this mechanism is available near Pinyon Rat (see lower 
portion of map in Figure 2.5) where Quaternary terrace deposits are isolated on older rocks. 3) 
Finally, the surfaces may have formed locally by erosion to a new base-level. Such a surface of 
low relief on exposed bedrock, when formed in a semi-arid region and situated between an allu- 
vial valley and a mountain, is commonly called a pediment (see, for example, Cooke and War- 
ren, 1973). The concave fo;ni of the topography (towards the SW), the thin veneer of weath- 
ered material, and the size of the area are all features in common with pediments. Absent 
from the Pinyon Rat area is a base-level defi led by an alluvial plain. 

This unusual situation may be explained by considering a mechanism presented by 
Wahrhaftig (1965). Wahrhaftig observed that exposed granites weather much more slowly than 
those rocks which are buried, and hence are moist for longer periods of time. Perhaps the rela- 
tively unweathered outcrops exposed on the southwest edge of the flat tend to act as a local 
base-level. Another possibility is that the change in lithology near the edge of the flat is 
responsible for reducing the erosion in that area. 

Whatever the cause or combination of nuses responsible for the creation of the surface, 
the most significant process currently affecting the surface layers is weathering. The near sur- 
face (at a depth of .5 m) is subject to a yearly temperature range of 18°C, with daily variations 
on the order of 3°C, On rare occasions ground frost will reach a few cm in depth. The aver- 
age annual rainfall has been 25.4 cm, for the period 1971 through 1980. Normally this 
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precipitalion cK'curs in ihe form of infrequent (but intense) thundershowers in the months ot 
August and September, along with less substantial prolonged winter storms. The top 1 m of 
ground at the site is nearly fully decomposed granodiorite. Below this level the weathered rcKk 
grades to highly competent grus at a depth of roughly 3 m. This material has decomposed in 
situ; dikes and other features of the rock are readily identifiable. At this depth standard soil 
moving machinery is no longer effective in excavating material. However, once a section of 
rock is removed it is found to be friable. The majority of the reference monuments at PFO 
were emplaced in the material at this depth. Conventional benchmark rods reach the depth of 
refusal in this layer (1-3 m). Evidence provided by drilling suggests that below 3 m to a depth 
of 25 m the material grades from grus to grus with corestones, and finally to jointed granodior- 
ite as illustrated in the right half of Figure 2.6. 

Wyatt (1982a) conducted a shallow seismic refraction survey at PFO which supports these 
observations. Figure 2.6 presents the possible range of the compressional seismic velocities for 
two locations at Pinyon Flat. These velocity corridors were determined using a procedure 
described by Garmany et J. (1979). The heavy solid lines are the extremal bounds on the 
solution to the seismic travel time inverse problem at PFO, based on the limited and inexact 
nature of the field observations. The dashed lines are the bounds at Edge. All possible models 
of the velocity profile lie within these bounds. The light solid lines indicate particular models 
that may be used to generate the observed travel times. It is unlikely that the layering shown 
in these mcxiels exists; rather the velocities probably increase smoothly with depth, as the 
effects of weathering diminish. 

There is a marked differenc’ between the seismic velocities at PFO and Edge. Velocities 
less than 2 km/sec, normally associated with sediments, are limited to depths of 8 m at Fxige, 
but may extend as deep as 18 m at PFO. The PFO profile is representative of the majority of 
the area, while the Edge data are more characteristic of the material at the limits of the flat, 
where outcrops (tors) are prevalent. Evidently the weathering is less near the edge of the flat, 
supporting Wahrhaftig's hypothesis for the formation of such surfaces. 

Below 20 m the velocities at both kx:ations tend to values that arc characteristic of grano- 
diorite subjected to low confining pressures ( .5 MPa) (Birch, 1966). Analysis of a core sample 
from PFO established the seismic velocity at a depth of 38 m to be 4.8 km/sec (C. Sondergeld. 
personal communication, 1980). This value is in excellent agreement with the mean velocity 
obtained by extrapolating the linear trend of the extremal bounds in Figure 2.6, and is 
representative of crustal rex'k velocities We may assume that instruments not adequately cou- 
pled to the r(Kk below the material with low velocities (~'20 m at PFO) will be sensitive to the 
displacements inherent in the weathering process. 

2.2 Instrumentation 

While PFO tixlay contains some 20 independent sensors, and experiments run by a 
number of different organizations (UCLA, UCSB, L-IXjO. NASA/JPL, AFGL, Leighton and 
.\ss(Kiates. CIRF^S. National Geodetic Survey. Cambridge, and Carnegie Institution) (see Fig- 
ure 2 7), we will confine our comments in this report to the results from the long laser strain- 
meters (LSM) These instruments were the first to be installed at the observatory and have 
undergone continual refinements in an effort to '.crease their accuracy and lower their noise 
levels 

Berger and Lovberg tl970) described the basic LSM design which was built and operated 
near the UCSD campus in San Uiego. California The optical configuration of the strainmeters 
is that of a Michelson Interferometer with a 732 m arm (see Figure 2.8). The changes in 
length of this arm are thus measured in terms of the wavelength of the laser light We use a 
Spectra Physics Model 119 for the light source as it provides a uninnxlal single frequency out- 
put. Fhe laser radiation is. of course, highly momvhromatic. but the center frequency is deter- 
mined by the mechanical properties of the laser resonator. Hence, the long term wavelength 
stability IS dependent upon this resonator's mechanical ability which is probably grxxl only to 
about a part in 10'’ over long periods To achieve a stability sufficient to measure earth strain' 
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Figure 2 1 Plan vi<*\k of existing instrumenlation at PFO 
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Figure 2 8 N-S ar»d E-W strainmeters 
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the wavelength of the laser radiation was locked to the length of a fused silica Fabry*Perot reso- 
nator, contained in a controlled environment. The length standard, then, for these experi- 
ments, is a 30 cm tube of fused quartz. Each of the three strainmeters is independently stabil- 
ized, but periodically the Fabry-Perot stabilizers are exchanged between the three strainmeters 
and the records examined for concurrent changes in the secular strain rates. The lack of any 
such changes gives as a high degree of confidence that, at least to the levels of strain rate we 
report herein, the length standards do not contribute to the secular strain observed. 

Following this basic design, with a few minor changes, three of these devices were 
in.stalled at PFO, one each year in 1971, 1972 and 1973. The instrumental statistics for these 
are given in Table 2.1. 

Both the NS and NW instruments are installed along concave topography and so it was 
necessary to rcsieer the beam part way along its length in order to keep it close to the ground. 
There are two prisms mounted in a housing at the center point (the low point) of the vacuum 
tube Each may be independently rotated so that the beam may be re-steered from this point. 
During construction the beam is steered to the center point, through the air, half the vacuum 
tube constructed and then, with a window at the steering prism and half the length evacuated, 
the beam is steered to the retro- re fleet or pier again through the air. 

The fust two strainmeters, the N-S and E-W, utilized the end pier design as shown in Fig- 
ure 2.8. Holes about i m in diameter were drilled into the surface to a depth where the drill rig 
(a caisson rig) could no longer penetrate. Rectangular columns of black granite, quarried to 
size, were then set into these holes with a non-expansive grout cementing the bottom meter or 
so to the hole walls. The remainder of the space between the column and the hole sides was 
filled with styrofoam balls. The column with its optics on top was covered with an insulated 
box and the whole apparatus enclosed in a temperature controlled building. 

The NW-SE strainmeter has a somewhat different end design. In an attempt to alleviate 
the effect of tilt of the end columns, it was decided to use roughly cubical end piers and bury 
the entire end structures in a vault. Figure 2,9 illustrates the NW strainmeter layout. 

During 1973, quartz tiltmeters (Wyatt, 1982a) were installed on each of the six strain- 
meter piers for reasons given later. 

The tilt t ansducers consist of three pieces of fused quartz that form a three-element 
differential tv /acitive transducer. Two of the elements are clamped together and plated with 
platinum elc.trodes to form the fixed plates of the capacitor. The third piece is plated on both 
sides and is suspended by two fine wires between the two fixed plates with a separation of .25 
mm on each side. The effective pendulum length is 6,75 cm with associated free period of 1/2 
sec and it is heavily overdamped. 

This whole assembly is clamped to a ground section on the side of the strainmeter granite 
piers The unit is hermetically sealed and temperature controlled. 

One of these instruments is mounted on each of the six strainmeter end piers. The 
instruments are readily calibrated by pushing on the piers and noting both the tiltmeter output 
and the strainmeter output, thus giving a lilt produced strain coefficient for each column. Typi- 
cal motion induced at the top of the pier, with a strenuous push is 6 x 10 <Yn with a 
corresponding till for a 3m column of 2x10*’ rad. 

After corrections to the strain records were made for the local tilt of the piers as measured 
by the pier tiltmeters, significant changes in the strain were .still observed, correlated with heavy 
precipitation, as we shall discuss in Section 2.3. These perturbations were especially obvious in 
the corrected records of the NW-SE strainmeter. The inability of the pier tiltmeters to correct 
for the liKal motions of the cubical piers of this sirainme’er led us to the design of a device to 
measure the displacement of the piers relative to points some 24 m below the ground surface 
An "optical anchor" (Wyatt <■/ a/,, 1982) was installed at the northwest end of the NW-SF 
strainmeter early in 1980 
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Table 2.1 

Strainmeter Parameters NS. EW ?.r.d NW 


Component 

NS 

EW 

NW 

Azimuth 

180.28'’ 

90.44'’ 

135.26'’ 

Dip 

o 

p 

-0.73'’ 

0.7" 

Length (m) 

732.2 

730.8 

730.4 

Height of end piers (m) 

3.8 (N) 

4.7 (W) 

1.2 (NW) 


2.9 (S) 

4.7 (E) 

1.7 (SE) 

Depth of end piers (m) 

1.2-2. 7 (N) 

1.0-3 2 (W) 

3.4 (NW) 


.8-1.8 (S) 

1.0- 2,7 (E) 

3.3 (SE 

Gain (strain/i.e.) 

1.7610 X 10 

1.7713 X 10-'" 

1.7666 X 10 

Filter time constant (s) 

507 

506 

521 



northwest - SOUTHEAST LASER STRAINMETER 
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Figure 2 9 NW - SE strainmeter 



22 


2.2.1 Optical Anchort 

Figure 2.10 is a schematic diagram of the optical anchor. A laser interferometer is used to 
monitor the difference in length between the two paths PqP\ and The horizontal shear 
strain can be computed from the change in this distance, provi.-^ed that the points P\ and Pi aie 
either well coupled (df| *• ? 2 ) or are undergoing much smaller displacements than Fq- Geologic 
evidence suggests that the near-surface at PFO is much less stable than the material at depth 
(Wyatt, 198.’a), so that both of these assumptions are likely to be true. 

Formally, we may describe the measurement in terms of a coordinate system in the 
bedrock below the instrumentation (Figure 2.10). The x axis is aligned in the direction of the 
shear measurement, the xi axis is in the vertical plane of Pq, P^ and /* 2 , end the x^ axis is per- 
pendicular to both. For this instrument, where the optical path lengths are much larger than 
the displacements to be monitored, the quantity recorded is: 

where and ?2 are the displacement vectors for the points Pqi Pi. and Pi respectively, 

and P| and Pj are the unit vectors along the two optical paths. To simplify this expression we 
assume that (do-Ji) ■ (do - d{i - d. Fnis is justified when the end points P\ and Pi are 
located in material more stable than Po (Id2l. Id 2 1 « Idol). Given that 02 " n - 9\ equation 
( 1 ) reduces to 

£) - 2 cos 0) X d • X| (2) 

Finally, for 0] ■■ 60®, equation (2) becomes 

Z) - d • X, (3) 

which is simply the horizontal displacement of point Pq relative to "fixed" points Pi and Pi. 

Figure 2.11 illustrates the arrangement of the optical anchor in the northwest vault of the 
Northwest-Southeast (NW-SE) strainmeter at PFO. Two boreholes, inclined at 30® from verti- 
cal, were drilled some 26 m from a point 2.5 m below the surface. The holes were drilled in 
the vertical plane aligned with the azimuth of the strainmeter axis, and intersected at a point 
roughly 1 m above the vault floor. A cement monument approximately 1 m^ in volume was 
formed at the top of the boreholes and a granite optics table placed on this pad. Reflectors 
mounted near the bottom of the boreholes served as the fixed reference points for the instru- 
ment. 

Proper orientation and alignment of the two boreholes were essential to the operation of 
the instrument. The driller (Rex Anderson Corp., Julian, California) was successful in the ini- 
tial attempt at this operation. A top-head drive drilling rig was positioned above the strainmeter 
vault, locked at an angle of 30® from vertical, and the holes were bored with an 16.5 cm air 
rotary percussion drill. To maintain a straight hole the hammer bit was followed by a 6.4 m sta- 
bilized drill rod (The stabilizer was made by welding three steel bars in a spiral along a stan- 
dard drill section so that the overall diameter was 16.5 cm.) The profiles of both holes had 
nearly uniform curvature, with a deviation of less than 5 cm from the straight line connecting 
the end points. The final clear aperture, necessary for the laser to illuminate the center of both 
remote reflectors, was 3,5 cm, barely larger than the beam diameter (2 cm). Table 2.2 presents 
other phy.sical and instrumental parameters. 

The most difficult part of the instrument design was devising a method to attach the 
remote reflectors to the sides of the boreholes in a stable manner. Although positioning an 
optical component with a stability of 40 ^m/year (corresponding to a strain rate of 50 
nstrain/year on the 731 m instruments at PFO) is easily accomplished in the laboratory, the 


♦This seciion is taken from Wyatt ei al (1982) 
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Fillure 2 10 Sthemaiic of the opiual unchor Points P^ and P-< arc locaterl at a depth of 24m P, And Ajare unit vec 
tors along the two optica! r«:hs (after Wyatt it «/. 1982) 
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f'lgure 2 1 1 Dei4il» o( the op»tc«l irnhor imerferonwier «rxj the reference momeni at ihe t’orih wesi end o( the 
Northwest Southwetl »ir»inmeler. Pihon Rti Obscrviiory (after Wyatt era/, 1912) 




Table 2.2 

Optical Anchor Physical and Instrumental Parameters 
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Physical Parameters 

Length of optical paths 
Path length difference 
Depth of interferometer 
Depth of retroreflectors 
Borehole angle relative to vertical 
Borehole azimuthal alignment error 
Laser beam intersection aperture 
Borehole inside diaameter 
Casing outside diameter 
PVC casing wall thickness 
Evacuated tube outside diameter 
Evacuated tube wall thickness 
Bellows outside diameter 
Bellows spring constant 
Bellows position imbalance 
Pressure on casing during cementing 


Instrumental Parameters 

Pressure in evacuated tube 
Pressure variations in tubing 
Laser output power 
Expanded laser beam diameter 
Beam divergence over 50m 
Temperature coefficient 
Maximum temperature fluctuation 
Precision of Optical Anchor 7. 

Cs'-bration accuracy 


25 m 
.5 cm 
2.44 m 
24.09 m 
30* ±2* 
±3* 

3.5 cm 
i6.5 cm 
14.4 cm 
7.1 mm 
10.2 cm 

1.6 mm 

5.7 cm 
lO*N/m 

3 mm 
400 kPa 


20 Pa 
5 Pa 
1 mW 
16 mm 
2.4 mm 
10 “ m/T 
2 “C 
9 X 10 * m 
.2 % 
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task is quilt difficult in the field. On the time scale of years, creep of the anchoring assembly 
in response to stress must be considered 

Figure 2 12 presents the detailed cross-section o^ the anchoring arrangement at the bot- 
tom of the inclined holes Each of the boreholes was cased with PVC pipe, with a 1.2 m stain- 
less steel anchoring assembly threaded onto the end of the casing. The assemblies were 
machined with many circumferential grooves and were cemented into place using non-shrink 
grout. Rubber pads were attached to the bottom of the anchors to reduce axial loading caused 
by borehole rebound and deformation of the cement. The grout was pumped through check 
valves in the bottom of the stainless steel assemblies and up the outside of the PVC casings 
until it reached the surface. During this operation the casings were filled with water to prevent 
the PVC pipe from collapsing. After the grout hardened, the water was evacuated from the 
casings by injecting compressed air. 

Near the upper end of the 1.2 m anchoring assembly, a tapered and threaded insert was 
machined to guide and secure the retro reflector housing. This housing was attached, by a pair 
of stainless steel bellows, to the end of a long stainless steel vacuum pipe. The bellows were 
required both to isulate mechanically the reflector housing from the vacuum pipe and to allow 
compensating ^lessure on both sides of the reflector. The vKuum pipe was supported only at 
the surface, assuring that the reflector housing was not subjected to exce^ive longitudinal 
forces (due to the weight of the tubing). Details of the retro reflector housing and the anchor- 
ing assembly are presented in Figure 2.13. 

Two detectors were monitored during the installation of the vacuum pipes and reflectors. 
One transducer, mounted on the reflector housing, indicated when the unit was fully screwed 
into the anchoring assembly The second detector was used to adjust the position of the 
vacuum pipe at the surface so that the forces of the bellows were balanced. 

A well pump was designed to keep the *nside of the borehole casings free of ground 
water, though this has proven to be unnecessary The pump consisted of an 8 mm copper tube 
that extended from the surface to a point below the reflector housing and a second tube (6 
mm) that connected to the first at a position 1 m above the bottom. The pump was operated 
by evacuating the longer tube, with the small diameter tube closed at the surface, so as to draw 
water up the larger tube to a height of roughly 10 m. Opening the second tube to the atmo- 
sphere forced most of the fluid above the joint to the surface. (Surface tension of the water 
kept the air from simply bubbling through the liquid.) Repeated "burpifig" of the system would 
lower the water level inside the casing. 

The differential path length measurement is easily accomplished using standard optical 
techniques (sec lor example. Strong, P58). Because the optical paths are nearly equal, the 
Michelson interferometer docs not need a laser with a very stable frequency. The frequency 
stability of an inexpensive Hc-Nc laser (A/// — 10 *) gives i maximum error of only 
5 X 10 '* m for the actual path length difference of .5 cm. The sense of the interference fringe 
motion was monitored by deliberately retarding, by 1/8 of a wavelength, opposing qucdranis of 
one of '.he optical paths, and detecting the resulting phase shifted fringe pattern (Berger and 
Lovberg. 1970) A reversible electronic counter was employed to record changes in length 
corresponding to 1/8 of a wavelength of He-Ne laser light (7.91 x 10 * m) 

To reduce possible errors, the pressure in the optical paths was lowered by con.nectmg 
both pipes to the vacuum pipe of the nearby laser strainmeier The pressure inside the stram- 
meter vacuum tubing was main.f^ined at I Pa (Berger and Lovberg, 1970). A rcstnction m the 
vacuum conduit between the strainmeier and the optical anchor assured stability of the strain- 
meter vacuum pressure. Tlie resulting pressure in the optica! anchor vacuum paths was approx- 
imately 10 Pa: more than adequate to eliminate temperature differentials in the two paths as a 
source of error However, pressure differentials jf only 1 Pa correspond to a displacement of 
8 X 10 * m. roughly the limit of resolution of the instrument. Maintaining the absolute pres- 
sure bdow 20 Pa limited the magnitude of the differential pressures to an acceptable level 
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Figure ’12 Anchoring assembly, shown cemenied near the bottom of one of the angled boreholes (after Wyatt et al., 
1982) 




DETAILS OF REFLECTOR HOUSING AND ANCHORING ASSEMBLY 



Figure 2 13 Details of the reflector housing and the anchoring assembly (after Wyatt et ai). 
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The completed instrument had a temperature coefficient of (approximately) 10 ** m/'C, 
probably from temperature sensitivity of the optical mounts on the granite optics table. No 
attempts were made to improve this, because the strainmeter vault was temperature regulated, 
ameliorating the problem. 

The optical anchor was calibrated by comparing the output of the system with the signal 
from the coaxial NW-SE laser strainmeter. Slight errors in the alignment of the optical paths 
were the main source of uncertainty in the gain of the instrument. Displacement of the strain- 
meter end monument provided an immediate calibration, as this displacement was measured 
(by the NW-SE strainmeter) in terms of the wavelength of He-Ne laser light, known to more 
than five significant figures. In practice, the accuracy of the calibration was limited by our 
reluctance to displace the reference monument so much that the bond with the earth might be 
broken. The final transducer gain was found to be 7.9 x 10^* m/count. 

2.3 Results 

In this section, we will report on the results of the laser strainmeters and their ancillary 
instrumentation. R ;sulis from the superconducting gravimeter are reported in Sections 4 and 5 
of this report. 

One of the firsi questions addressed in the analysis of the PFO strain observations, and a 
general question regarding all strain observations, is how do we know that observations on the 
surface at one place are representative of the regional crustal strain? 

Distortions to the regional strain field may arise from inhomegeneities in material proper- 
ties in the neighborhood of the observing point. Topographic irregularities will also cause an 
effect, distorting the strain field from that which would be observed on a smooth surface. 
Berger and Beaumont (1976) investigated the distortions lO a homogeneous strain field caused 
by the local topography at PFO and a crude two-dimensional model of the regional geology. 
This study concluded that a 10-20 percent variation in the measured linear extensions might be 
expected, compared to those on the surface of a homogeneous ellipsoidal earth. Clearly, as we 
move from place to place, even if the applied strain field were homogeneous, we would expect 
fluctuations on (hat order. 

To test if our point observations of strain come close to matching our theoretical models, 
we can use the earth tides to "calibrate" the observatory. Because of the well understood theory 
of the tidal potential and a quit? accurate spherically averaged earth structure, it is straightfor- 
ward to calculate the earth tide strains or tilt at any location for any azimuth. In the real world, 
however, perturbations to this theoretical earth tide are caused by ocean loads applying an addi- 
tion signal at the same frequency but not in phase with the potential as well as the distortions to 
the strain field due to local inhomogeneities such as topography and geology. 

In Section 5 we review the work of Berger and Wyatt (1973), Berger and Beaumont 
(1976), and Agnew (1979), who analyzed the .strain tides observations at PFO and made calcu- 
lations of the ocean load perturbations. Figure 2.17 shows the results for the semidiurnal tidal 
lines (.Agnew. 1979), While there a'e certainly some errors in the models of the ocean tides 
used to calculate the ocean loads, it is likely that even larger errors are caused by inadequacies 
in the modeling of the local inhomogeneities. Nonetheless, the analysis gives us confidence 
that the observed strains are within 10-20 percent of our theoretical values and we have no rea- 
son to suppose that this conclusion cannot be extended to other frequencies. 

There is ano'her naturally applied strain field that we believe we can model fairly accu- 
rately and thus use as a test of our observations. This is the strain field that is produced when 
an earthquake occurs We are considering here the static dislocation of ihe fault movement 
rather than the dynamics of the elastic wave generation. 

Wyatt (1982b) examined the tilt and strain records recorded at PFO for the October 25. 
1979. Imperial Valley earthquake. Although the theoretical displacements for the earthquake 
utilizing disUxiation models are clearly not as accurate as the tidal modeling, the observations 
agreed to within 25 percent of the model with basically one free parameter, the seismic 
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moment, ailjusied to fit three independent observations. Figure 2,14 shows the offsets recorded 
on two strainmeters (the third meter was rendered inoperative by the seismic accelerations 
which misaligned the laser beam). Figure 2.15 shows the residuals after removing the tides and 
the strains due to the dislocation model of the earthquake. 

Both the tidal and seismic dislocation results give us confidence that, at least at periods of 
a day and less, the strain field we are observing at PFO is a quite accurate measure of the 
regional strain field. 

As we extend our measurements to longer periods, however, and attempt to observe the 
secular strain over months and years we must face the questions of instrumental instability and 
local displacements that will cause noise in the records. 

The initial design of the laser strainmeter to which the North-South and East-West com- 
ponents were built used long (—5 m) columns of block granite as fiducial points or end piers. 
These were set into holes drilled some 3 m into the surface material as described in Section 2.2 
(see Figure 2,8). It was hoped that by this construction, the instruments would be detached 
from the surface layers that might be subject to environmentally caused displacements. This 
proved to be a false hope. Whenever there was a rainfall of over 2.5 cm. there would be large 
apparent strains recorded on both of these meters. A clear correlation with air or ground tem- 
perature was not apparent. 

In an attempt to correct for these environmentally induced pier displacements, tiltmeters 
were attached to the piers at either end and the tilt along the axis of the strainmeter monitored 
(Wyatt. 1982a). Figure 2.1o shows the results for a six year section of the East- West strain- 
meter The top curve shows the observed uncorrected strain. The next two curves are the pivr 
tilts (converted to strain units). The pier tiltmeters are calibrated in a manner which allows for 
direct removal of the displacements due to tilts about the base of the end piers. The curve 
marked "Corrected Strain" in the figure is obtained by using this calibration when applying the 
Fjtst and West pier tilt corrections to the raw data. The lower curve demonstrates the improve- 
ment in the results when the pier tilt records are least squares fit to the strain record. The vari- 
ance in the residual is decreased by nearly a factor of four and the value of the free parameter 
m the fit implies that the piers are rotating at a point significantly below their base. 

The third of the laser strainmeters to be built was designed differently. As described in 
Section 2.2. its end piers consisted of nearly cubical granite blocks set on the floor of under- 
ground vaults. The base of these piers was approximately 3.5 m below the surtace. The idea 
behind this design was that by eliminating the long thin columns, the sensitivity to pier tilts 
should be greatly reduced. The observed uncorrected records showed just this. Unfortunately, 
after the corrections for pier tilts were applied to the NS and EW cornponenta. the residual 
strain was apparently more free from environmental artifacts than the NW-SE. On that com- 
ponent, the pier lilts, not unnaturally, made little difference but yet the corrected strains 
showed an unrealistic secular trend in addition to usual rainfall effect. 

As we discussed briefly in Section 2.2, the inability of the pier tiltmeters to improve the 
NW-SK strain records led us to install the optical anchor on the NW pier, late in 1979 Shortly 
after its installation, a Mi - 5 3 earthquake occurred on the San Jacinto fault at an epicentral 
distance of only 14 km from the observatory. Figure 2.17 shows a portion of the strain records 
for the NS and EW instruments around this time period. The offsets in the strain records have 
been removed but the residuals clearly show a strain "precursory" to the earthquake which we 
belic\e was caused by the heavy precipitation that occurred. (Just prior to the earthquake. 22 5 
cm of ram fell on PFO in a nine day interval. The average yearly rainfall over the past nine 
years has been 25.4 cm ) Figure 2.18 shows the NW-SE strain both observed and corrected by a 
combination of the NW optical anchor and the SE pier tilt. The optical anchor record reveals 
that almost all of this anomaly was caused by a local displacement of the NW pier. The residual 
after removal of the optical anchor signal is remarkably smooth so we are able to see fluciua- 
i.ons on the order of 10 strain/day! Figure 2,19 .shows the residual strains on all three com- 
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F igure 2 ! “i Possible post scismit deformation recorded by tFie North-South and Northwest-Southeast strainmeters at 
PFO The lower traces were produced by removing the tidal components from the strain observations (after 
Wyal! flat. 19821 
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Figure 2 l^ Records frorn ihe 7Jlm East-West sirainmeter at PFO, indtttMiii ihe pier lilt vsries used to Cf.ect the 
stram obscrsations (converted into equivalent displacements), and the sirainmeter dis|iiaccmcnis corrected in two 
different m.tnners Mso presented is the record of the daily preopitation 
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Kigure 2 17 S«r»in records for NS and EW instruments before the - 5 3 Horse Canyon earthquake 
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NORTHWEST - SOUTHEAST STRAI N METER 

PFO 




Figure 2 18 NW-SE strain observed and corrected by a combination of the NW optical anchor and SE pier tilt before 
the W/ - ^ 3 Horse C'anvon earthquake 
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poncnts for comparison. Wc anticipate that the corrected NW-SL strain signals will be reduced 
below the level of the current optical anchor corrections ( — 10 ’ strain/year) once the SE opti- 
cal anchor is installed and operating. 

Figure 2.20 shows the results of applying the optical anchor corrections over the entire 
period of operation. To understand fully this data, and (he longer record of (he NW-SE strain 
presented in Figure 2.21, we must explain (he following; 

1 When the optical anchor was installed, the original granite strainmeter pier at the 
NW end was removed and a cement monument replaced it. (As the strain measure- 
ments at this end arc now referenced to the points at a 24 m depth via the optical 
anchor, the surface pier need not be dimensionally stable.) We now believe that the 
secular strain rate up to this time was caused by a slow movement of the granite pier 
(which was not well constrained at its base) downslope towards the SE end. 

2. Having discovered this, additional cement was added around the base of the SE pier 
at the start of 1980. This is reflected by the sharp inflection in the strain record at 
this time. Following this alteration, the SE pier's immediate response to rainfall was 
considerably reduced so that very little effect of the heavy rain in early 1980 can >e 
seen in the strain record 

From 1980 onwards, the principal feature of the strain record is the bump near 1980.5 and 
the rise during late 1980. Neither of these two features is removed by the optical anchor 
correction and represents either continued SE pier motion (soon to be measured by the new SE 
optical anchor) or a real strain. Figure 2.21 shows the corrected strain observed at PFO over 
the past seven years. Also plotted in this figure are the strains deduced by analysis of the 
USGS Pinyon Trilateration Network (Savage, personal communication, 1981) and the local 
rainfall. 

Perhaps the most striking feature of the PFO strains is the large NW-SE secular compres- 
sion from 1974 to 1978. However, as we have explained in Section 2.2, we believe this result 
is wholly instrumental, being caused by a local movement of the end piers. The data on this 
component should have improved after the installation of the optical anchor on the NW end in 
late 1979 (The optical anchor on the SE pier is under construction at the time of this writing. 
October, 1981.) 

The secular change in the NS and EW components is much less. The agreement with the 
geodetically determined strains is remarkably good at least up until late 1978. At that time, 
there seems to have been a change in the normal strain pattern reflected in both sets of data. 
Indeed, a change in the strain pattern in all of southern California was noted to occur at about 
this time (Savage et al„ 1981). It is not at all unlikely that (he differences between the geodcti- 
cally determined strains and those from the strainmeters of PFO represent spatial inhomo- 
geneities in the strain held. 

Figure 2.22 shows the results of the USGS trilateration network called ANZA and a sub- 
network referred to as Pinyon (see Figure 2.4). Tnc former is a point network spanning the 
Elsinore. San Jacinto, and San Andreas faults while the latter is an 1 1 station network centered 
on PFO. In theory, if the assumptions of homogeneity and temporal stability that are used to 
reduce the line length changes (Sevage ef al„ 1981) to strains are valid, then one would expect 
these networks to give the same strains with the smaller network having larger error bars. 

In general, the results of the two networks do agree but the differences between the two 
may reflect spatial inhomogeneities in the strain field on the order of the observed strains 
Thus, we are not surpri-^d to sec differences in the strains determined at PFO and those of the 
encompassing trilateration network. 

.Mmost all the major fluctuations of the PFO strain records about their linear trends can 
be attributed to the local effects of rainfall. There may be one or two fluctuations about the 
trends not coincident with a local rainfall that can see be seen upon close examination of the 
records but these are of small amplitude and dubious significance We cannot explain the 
apparent change in response to rainfall of both the North-South and Ea.st-West components 
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occurring in late 1978. It seems unlikely that this is connected with the regional change in the 
strain held mentioned above but may very well be related to the change both locally and in all 
of southern California from an epoch of drought to one of more normal precipitation. 

In summary, using the deformation data from PFO observed over the past seven years we 
can; 

1. Compare deformation observed at PFO with deformation deduced from a trilatera- 
tion network which encompasses the observatory. 

2. Determine upper limits on the rate of change of deformation (strain and tilt) to be 
expected. 

1 Comparison of strain. As we have shown from the strain observations at PFO over 
tf . period 1974-1981 and the results from the USGS Pinyon Trilateration Network, the long 
" sect lar strain recorded in the North-South and East-West directions agree with the trila- 
tu.iUon ' suits as well as might be expected considering the difference of the two techniques. 
The results from PFO are within two standard deviations of the trilateration net data. The 
Northwest-Southeast strainmeter, however, has not produced consistent secular strain rates. 
We believe the reason for this has to do with the details of the attachment of the instrument to 
the earth (Wyatt, et al, 1981). Recent additions to this instrument have improved the charac- 
teristics of the records greatly. 

2. Rate of strain. With regard to the time variations of the strain, almost all of the obvi- 
ous fluctuations in the three strain records are correlated with local rainfalls of 2.5 crr= and 
larger. Wyatt (1982a) has made a thorough study of near surface displacer, ents and concluded 
that these signals are artifacts of the local motion of the end monuments-the strainmeter 
attachments points. Thus, these records must be taken as o 'erestimates of the normal back- 
ground strain noise. Perhaps the most succinct way to summarize this data is to calculate the 
averaged power spectrum of the three records. Figure 2.23 shows the results. Table 2.3 gives 
the average standard deviation in a one decade bandwidth centered at a set of discrete frequen- 
cies. (This bandwidth was chosen to obtain agreement with intuitive estimates of the "ampli- 
tude" spectrum for a power spectrum whose frequency dependence is approximately 1/ /^.) In 
general, if we eliminate artifacts obviously caused by rain, seismic and tidal signals, we have 
seen no strain changes at periods shorter than the average repeat time of the trilateration nets 
greater than 7 x 10“^ that appear on two or more instruments. (While it is possible to produce 
a deformation that changes only one component of strain, it is extremely unlikely to occur in 
nature. Hence, we label as instrumental or environmental artifact those strain changes seen on 
one instrument only.) 
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Table 2.3 

Average Standard Deviation in a One Decade Bandwidth centered 
at a Sc i of Discrete Frequencies 


Frequency Power Spectral 
Cycles Density 

Per Day Relative to eVHz 


Standard Deviatiion 
In a 

One Decade Bandwidth 


001 

- 68 dB 

1.8 X 10 ' 

.01 

- 85 dB 

2.8 X 10- * 

.1 

-112dB 

3.3 X 10 ’ 

1.0 

-140 dB 

5.4 X 10 

10.0 

-162 dB 

1.4 X 
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3.0 S«Uinicity Studies Related to the Gulf of California 

During ihis projeci we have had close ties with Mexico in connection with operation of 
the Northern Biya California seismic array and other seismicity studies in Mexico. One Mexi- 
can graduate student in seismology, Luis Munguia, is currently completing his Ph D. Alejandro 
Nava completed his Ph D. in March, 1980. Alfonso Reyes, Chief of the Geophysics Depart- 
ment of CiCESE (Centro de Investigaciones Cientifkas y de Educacidn Superior de Ensenada), 
did his graduate work at IGPP and is now completing his Ph.D. thesis. 

We have recently completed studies of the seismicity of Northern Btya California (Brune 
ft al„ 1979). the structure of Northern Bcya California (Nava and Brune, 1981a); the source 
mechanism of two Northern B^a California earthquakes (Nava and Brune, 1981b); aftershocks 
of the Oaxaca earthquake of November 29, 1978 (Reyes ei ai. 1980); spectra of aftershocks of 
the Oaxaca earthquake (Munguia et ai, 1980); and the location of the hypocenter of the 
Imperial Valley earthquake of October 15, 1979 (Chavez ei al., 1981). At the present time, we 
are preparing papers on a recent swarm of earthquakes which occurred in the Laguna Salada 
region (Reichle and Simons, in preparation, 1981), on aftershocks of the Petatlan earthquake of 
March 1979 (Suarez and Reichle, 1979), and on the source mechanism and aftershocks of the 
recent Victoria. Baja California, Mexico 'u 'hquake of June 1980. 

Brune vi al. (1979) summanzed the current evidence relating to seismicity and faulting in 
Northern Baja California. The area is one of extremely high seismicity, at least 13 earthquakes 
of magnitude greater than six have been recorded since 1900. Activity on the San Miguel Fault 
IS particularly high. Relocation of the 1954 Bi^a California earthquakes indicates that they were 
on the San Miguel Fault (and not on the Agua Blanca Fault as had been suggested earlier). On 
.August 19, 1 978, a small earthquake (M— 3.5) occurred on the San Jacinto Fault zone just 
southeast of Tijuana, emphasizing the possible hazard to Tijuana and San Diego. The region of 
the Mexicali V'alley is the locale of the main plate boundary and the seismicity there is the 
highest in the region, posing a severe seismic hazard to the city of Mexicali. 

Recently, Nava and Brune (1981a) completed an explosion-earthquake reversed refraction 
profile in the Northern Peninsular Ranges between Corona. California and Valle de Trinidad, 
southeast of Ensenada. They found a relatively high velocity crust (6.6 to 7.0km/sec) and a 
crustal thickness of ~40km. 

Nava and Brune (1981b) also compared two earthquakes, one from the Salton Trough 
region and the other from the granite Peninsular Ranges. They found that the main difference 
boiween the two earthquakes was the complexity and much higher relative excitation of surface 
waves by the Salton Trough event as compared to the relatively simple nature and low excita- 
tion of surface waves for the Penninsular Ranges event 

In a study of the Oaxaca earthquake of November 29, 1978, Reyes ei al. (1980) analyzed 
the aftershiKk sequence distribution along with leleseLsmic data. The overall analysis .suggests 
that the rupture process a.ssociated with this earthquake cannot be interpreted in terms of sim- 
ple plate tectonic processes. One interpretation explains the fact that the heaviest damage was 
m the high elevation regions; the rupture probably started at the intersection of two faults, pro- 
pagating toward the trench axis and therefore, focusing its energy along this direction. The 
increased stress field due to the propagation of the initial break probably increa.sed the stress on 
the secondary fault triggering a second event that most likely focused its energy toward the high 
elevation regions. This secondary fault probably stopped the main fault from extending into 
the deepest regions of the nijyor plate boundary. However, this explanation does not account 
for the lack of damage in the lowlands along the coast, where the heaviest damage would be 
ox[iected because o( amplification effects in the low rigidity sediments A second interpretation 
consistent with the teleseismic data was formulated. The main rupture started along the secon- 
dary fault propagating parallel to the trench axis from the deepest part of the fault toward the 
high elevation regions This propagation obtained from the P wave teleseismic data, but it adds 
complexities to the rupture processes along the subduction zones. 
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Mungui'a ei al. (1980) are currently studying a swarm which occurred in March, 1978 near 
the town of Victoria, Northern Baia California. They found that many of the larger quakes 
(Kcurrcd al depths of 15km. Hypocenters in the Salton Trough, where the faults are usually 
covered by —5km of sediments, are usually relatively shallow (~10km) and in the large part 
confined to the sedimentary section. Thus, the depth of these earthquakes was surprising, as 
wee the apparent high stress drops for these events, some as high as 1 kilobar (typical stress 
drops for large earthquakes are —30-100 oars). We are currently studying this swarm and the 
study will constitute part of the Ph.D. thesis of Luis Mungui'a. 

Chavez et al. (1981) studied the October 15, 1979 Imperial (Mexicali) Valley earthquake. 
Ml - 6.6, which occurred in the vicinity of several permanent seismic networks operating on 
both sides of the international border. Because the event lay outside the individual arrays, ini- 
tial Ux'ations calculated using the data from any one group of stations were misleading, 
incorrectly placing the epicenter in the United States. By pooling the arrival times from all sys- 
tems. however, complete azimuthal coverage was obtained and a reliable hypen 'ntral solution 
was produced. The mainshock epicenter was located 3km south of the international border, 
approximately 10km east of Mexicali, B^a Calilornia Norte, Mexico. The 9.96kni depth 
obtained is slightly greater than the 6 to 8km average depth for Imperial Valley earthquakes 
(Johnson, 1979). The significance, if any, of this is not well understood at this time. The 
greater depth can help explain, however, the fact that although the earthquake took place 
within Mexico, the only surface rupture observed was north of the border. 

Reichle and Simons are preparing a paper for publication on a lengthy earthquake swarm 
(December 1975 - June 1976. several M— 5) which occurred near Laguna Salada, Northern 
Biijd California. Mexico. The swarm is being studied with California Institute of Technology 
and the then existing Northern Baja California station data. Since this swarm did not occur 
along the most active branches of the Plate boundary (the Imperial-Cerro Prieto Faults), it may 
provide insight into the transfer of seismicity from the Gulf of California system to the various 
Southern California faults. 

t')n 9 June 1980, at 03:28:20 (GMT), an earthquake of magnitude 6.2 (.5// ) occurred in 
Mexicali Valley, Northern Baja California. Mexico. The instrumental epicenter was roughly 
10km southeast of the town of Guadalupe Victoria, Four hours after the main event, the first 
of a portable array of seven analog and 12 digital seismic event recorders, both from CICESE 
and IGPP (UCSD). began recording the aftershock activity. 

The last recorder was removed from the region on 5 August 1980. following the accumu- 
lation of over 50 smoked-pafier records and 60 digital cassettes containing aftershock data. The 
contents of all cassettes have been transferred to computer-compatible tape and time correc- 
tions have been determined for all stations. All usable signals have been edited from the 
computer-compatible tat>cs and archived. 

Investigators at CICESE have located many of the aftershcxks of the Victoria earthquake. 
Thev are clustered principally on a segment of the Cerro Prieto Fault close to the two commun- 
ities. \ ictoria and Pescaderos, which experienced the greatest damage from the mainshix'k. A 
number of events are also found scattered throughout the area between the southern end of the 
Imperial Fault and northern end of the Cerro Prieto Fault. The instrumenially determined epi- 
centes of the mainshrxk (based on initial P arrivals at nearby USGS and Mexican stations) is 
near the Cerro Prieto Fault, but is approximately 15km southeast of the region of aftershixks 
and major damage. A similar phenomenon characterized the 15 tVtober 1979 event on the 
Imperial Fault (Chavez er al.. 1981) 

Additional work at CICESE on the aftershocks of this event includes fcx'al mechanism 
studies and a comparison of this event with previous events on the same fault system. At 
IGPP. S-wave spectra are being computed from some of the largest aftershocks 
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An imix)riani aspect of the Victoria mainshock is the strong motion record obtained at a 
station in the town of Victoria. Vertical accelerations there exceeded Ig in both the positive and 
negative directions for several cycles, and the horizontal acceleration reached .85g. 

Soon after the Victoria earthquake, the region was inspected carefully from the air and on 
the ground by observers from both IGPP and CICESE. Many fissures and sand blows and 
other surface manifestations were seen striking roughly perpendicular to the strike of the Cerro 
Prieto fault system. However, no surface features paralleling the fault system were detected. 
This contrasts sharply with observations of surface faulting on the Imperial fault following the 
15 October 1979 Mexicali event. 
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4.0 Development of the Superconducting Gravimeter 

This work was motivated by data obtained with our superconducting gravimeter at PFO. 
That data indicated that variations in gravity of about ±5 n^tals could be observed over arbi- 
trarily long periods of time. Using funds from this grant, and several others, five instruments 
were built in 1978 for deployment to new sites. The object was to measure non-periodic varia- 
tions in gravity at various locations to attempt to measure vertical crustal motion and to learn 
what the limitations on such measurements might be. Ultimately, the results were to be inter- 
preted, along with spaceborne measurements of crustal motion, to give complete information 
on vertical and horizontal motion and associated mass displacement beneath the surface. 

Subsequently, instruments were operated at Otay Mountain (San Diego), Lytle Creek 
(north of San Bernardino), The Geysers (near Santa Rosa), Boulder, Colorado, and Greenbelt. 
Maryland. At Greenbelt. very severe problems with electrical interference were encountered so 
that only very short segments of useful data were obtained. The results from the other stations 
are summarized in Table 4.1. Plots of data after removal of tides and barometric pressure 
eft'ects are shown in Figures 4. 1—4.5. Some comments on the results are included in Table 4.1. 
We conclude that some "events" were observed on the gravity field at The Geysers, and a slow 
drift in gravity was observed at Lytle Creek until about two months before the instrument was 
turned off. The other stations reveal no significant variations in gravity within limits of ±5 
jugal except at Otay Mountain where large tilts produced artificial variations of ±10/xgal. 

At the time of these measurements all of the instruments were operated singly at isolated 
kxations so that it was not possible to prove that the different results were a consequence of 
different operating characteristics of the various instruments Our most recent work has con- 
centrated on this question by returning all of the instruments to the lab for simultaneous side- 
by-side testing. Some sources of spurious noise have been identified and eliminated and others 
arc still being diagnosed. However, no long term drift has yet been found. It is found that two 
properly operating instruments produce signals which are identical to within less than 0.2 /xgal 
This means that the ±5 ngi\ fluctuations are probably real changes in gravity which must be 
understood. It also means that long term observations, which average over these fluctuations 
(characteristic time of order one week), will determine steady gravity variations to better than 5 
Mgal. 

In summary, the project successfully measured gravity variations at several locations to 
within 5 /ugal. Subsequent work has shown that the superconducting instruments can measure 
gr; vny variations to within about ±0.2 Mipl 
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correction for tilt as determined by periodic relevelings of the instrument. The top trace 
the gravity residual after application of this correction. Exclusive of the winter months, 
where the tilts are greatest, the signal shows no raonotonic drift. 
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A ^79-day record from Eioulder, Colorado. The top trace is gravity after removal of the tides, 
the second is the barometric pressure, and the third is the residual gravity after removal of 
the pressure effect. 
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S.O Earth Tide Studies 

The time varying changes in the gravitational potential of the earth caused by the motions 
of the earth-moon-sun system give rise to the earth and ocean tides The earth tides are in 
equilibrium (or nearly so) with the tide potential as the highest frequencies in that potential are 
far below the lowest elastic eigenfrequency of the solid earth. The laterally averaged elastic 
properties of the earth are well kno«vn from seismological studies and hence the earth tides on 
such an occanless earth model may be readily calculated. The ocean tides, in contrast, are 
nonequilibrium tides as the ocean basins resonate near tidal frequencies. These ocean tides 
lo&d the earth's surface and cause significant perturbations in the observed earth tides. If the 
ocean tides were well known, it would be an easy matter to calculate these perturbations as the 
response of a laterally averaged earth model to such surface loads is known (Farrell, 1972a). 
Unfortunately, such is not the case. The global ocean tides are not well known nor easy to 
measure except along coasts. 

A second class of perturbations on theoretical earth tines arises from the departure of the 
earth from radial symmetry. The effects of geology (i.e. changes in rock type and hence elastic 
parameters), topography, and underground vaults if the instruments are placed therein, are all 
observable in the analysis of the tides of the real earth. Further, if temporal changes in elastic 
characteristics take place, as bas been suggested to occur prior to earthquakes, then the pertur- 
bations on the earth tides due to these changes may be observed. 

Our s'.udies of the earth tides then are centered in three areas; 

1 The ocean tides-improving and constraining global ocean tide models. 

2 Developing three-dimensional models of the elastic structure around tidal observa 
lories 

3. Investigating the temporal stability of nearby elastic structure and its connection with 
'he seismotectonics of the area. 

The response of the three basic kinds of instruments used to study the earth tides can be 
characterized as follows; 

1 Gravimeter or vertical accelerometer. This instrument responds to both vertical 
acceleration and, due to the earth's vertical gravity gradient, to vertical displacement. The tidal 
vertical acceleration is predominantly due to the direct lunar attraction On a perfectly rigid 
earth, this is all that would be observed. On an elastic earth, however, there is an additional 
effect that amounts to some 16 percent of the direct attraction and is additive. Thus the "gra- 
vimetric factor" is 1 16. 

The response of a gravimeter to surface mass loads likewise arises from the direct 
Newtonian atu-ction of the mass and from the vertical displacement at the observing point 
caused by the surface load. Figure 5.1 (after Farrell. 1973) illustrates the relative importance of 
these two effects as a function of di.stancc from the load. TTicy are equal at about 3° (330 km) 
and at greater distances, the direct attraction dominates, falling off like r '. 

2 The iiltmcter or horizontal accelerometer This instrument responds to both horizon- 
tal acceleration and tilt relative to the local vertical. In the case of :he earth tides these effects 
arc in opposition with the till (a result of the earth’s elasticity) being about 30 percent of the 
direct effect Thus the diminishing "tilt factor" is about 0 7. The response to surface mass 
loads IS dominated near the source by the elastic or surface tilt part of the response. At about 
5'" (.SOO km) It equals the direct Newtonian attraction However, the ratio of the horizontal 
Newtonian attraction to the vertical Newtonian attraction vanes like 

~ - -tan (e/2) 

s * 

so that at. say, 5' .V - .04, V' and thus the tilt is about 25 times the gravity perturbation. 
Close U) the source, the ratio of total gravity perturbation to total tilt is even smaller and in fact 
the tilt perturbation is almost solely due to the slope of the surface 
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Figure 5 1 The elastic and Newtonian parts of the gravitational potential Green function. The ampliluoes have been 
multiplied by (ofl/jtf) x !0'^ (a is 6 371 x l0^n. 0 is distance in radians, g is 9 8m s”^) to give relf ve ampli- 
tudes 
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3. The strainmeter. This instrument is a linear extensometer, measuring changes in 
length, AZ, along some base length L. Its response to the tidal potential is controlled solely by 
the elastic behavior of the earth for if the earth were perfectly rigid there would be no change 
in length. Lacking the direct effect, the strain tides are thus an order of magnitude smaller 
than either the tilt or gravity tides. This lack of direct response to gravitational potential attrac- 
tion may be deemed an advantage if one wishes to study the elastic structure of the earth, as its 
response to surface mass loads is also purely elastic. 

5.1 The Ocean Tides 

We believe that precision measurements of the body tides made with calibrated low noise 
instruments such as the super-conducting gravimeter and the laser strainmeters can significantly 
contribute to the modeling of the global ocean tides. Both these instruments are weil cali- 
brated, the gravimeter to better than 1/2 percent, and the strainmeters by their nature as fringe 
counting interferometers. It has been demonstrated that they both are limited by only earth 
noise in the tidal band (Warburton et al., 1975, Berger and Levine, 1974). 

As a first effort along these lines we have calculated the theoretical gravitational load tides 
of two major tidal constituents, M 2 (principal lunar semi-diurnal) and Oj (principal lunar diur- 
nal) for comparison with the observed tidal gravity signal at La Jolla and Pinyon Flat. These 
constituents were chosen because they have the best known global distributions, they are well 
removed in frequency from the effects of radiational tides, and there are no predicted anomalies 
in the earth tide admittance due to core resonance at these frequencies. 

The load tides were calculated by convolving Farrell’s Gutenburg-Bullen A gravity load 
Green’s functions (Farrell 1972) with theoretical models of the M 2 and 0\ global ocean tide 
distribution taken from Hendershott and Munk (1970) and Tiron et al. (1967) respectively. 
These global tide models are numerical solutions of the Laplace Tidal Equations (LTE) using 
coastal observations as boundary conditions. The models are in general agreement with the few 
existing deep ocean pressure observations (see for example Munk et al. 1970) and island obser- 
vations. However, the predicted tides in some areas ire at variance with the known ocean tide 
distributions, mainly because significant features in ihe ocean bottom topography cannot be 
described in sufficient detail on the coarse finite difference grids ( 4 "’ to 6 °) used in the calcula- 
tions. Smaller errors also arise because the theoretical results represent incomplete solutions of 
the LTE, and because the gravitational self- attraction of the tidal water mass has been ignored. 
More recent attempts (Farrell 1973, Hendershott 1972) to include these factors in the calcula- 
tions are inconclusive due to a numerical error. 

Fortunately, the ocean load tide in southern California is dominated by the tide in the 
northeast Pacific Ocean In this area regional mrvdels for the M 2 and O] tides based on coastal 
and ocean bottom pressure gauge observations have been developed by Munk et al. (1970). 
Their results were used to modify the M 2 and 0\ global tide models in the north Pacific Ocean. 
The modifications to the 0\ model are minor because the Tiron, et al. chart is in good (perhaps 
fortuitous) agreement with the observations. Figure 5.2b The modification of Hendershott’s 
M 2 cotidal chart was, however, more drastic and involved moving his amphidrome, located 
northwest of Hawaii, to the location proposed by Munk, et al. 1500 Km southwest of San 
Diego, Figure 5.2a. 

Calculation of the global load vector is complicated, however, 'xicause the Tiron, et al. 
and Hendershott models do not conserve tidal mass. Conservation of mass is unimportant in 
these solutions of the LTE where a semi -permeable coastline has been used to simulate fric- 
tional energy dissipation on the continental shelves. However, for the gravity convolution, tidal 
mass must be conserved (Farrell, 1972b). Because the distribution of non-conserved tidal mass 
is not known, an artifice such as distributing this mass equally over all the oceans must be used. 
This mass is small (7.01 gm/cm^ - 90.45® for M 2 and 0.47 gm/cm^ - 28.30® for 0|) but not 
negligible. In order to estimate errors incurred by this enforced mass conservation, results are 
presented for three different methods. In the first method the non-conserved mass is distri- 
buted uniformly over all oceans. In the second method each of the major oceans is forced to 
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conserve its own tidal mass. The two methods subtract significantly different tidal masses from 
some oceans, in particular the Indian Ocean, but the total load vectors are similar. These vec- 
tors should be compared with the results of the third method where no attempt is made to con- 
serve mass. We have no way of determining an accurate method of mass conservation but we 
assume that the errors incurred by this lack of knowledge are no larger than the ditference in 
the vectors computed by methods 1 and 2. 

The observational results presented here were calculated from three sets of gravimeter 
data. Runs 1 and 2 were taken at La Jolla from a 54 day and 100 day record. Run 3 was taken 
at PFO for a 365 day period. 

The complete data records were cosine tapered and Fourier analysed to determine the 
observed amplitudes and phases of the 0\ and M 2 tidal constituents. An identical analysis was 
performed on the theoretical rigid-earth tide which was generated for the same time period as 
the observed data. The results of the Fourier analysis are shown in Figure 5.3. The amplitude 
of gravimetric factor 16 1 is the ratio of the observed tidal amplitude to the theoretical amplitude 
and the local epoch is the phase difference between the observed and theoretical tides. The 
normalized amplitude I AS I and the local epoch of the "measured load vector” were calculated as 
the vector difference between the observed gravimetric factor 8 and an assumed theoretical gra- 
vimetric factor of 1.160 with zero phase lag (at the local equilibrium potential). The amplitude 
of the "measured load vector" is also given in fx gal and is calculated by scaling |AS| using the 
local amplitude of the theoretical gravity tide of the rigid-earth. 

The observational results can be compared with the theoretical load vectors and gravime- 
ter factors computed in a number of ways. The sum of the oceanload vector and the computed 
earth tide, described above, should equal the measured earth tide if the ocean load model is 
correct, and no other effects contribute to the loading. Using the results of method 1 (conserv- 
ing total mass of all oceans) the observed 0\ tide is larger than this theoretical tide by 0.15 per- 
cent at PFO and smaller by 0. 10 percent at La Jolla. Using method 2 (mass of oceans con- 
served individually) for Mi the observed tide is smaller by 0.20 percent at PFO and by 0.05 
percent at La Jolla. These differences are within our estimated calibration error. The 
differences between 0[ and M 2 , however, cannot be a consequence of calibration error. The 
calculation by method 2 appears better than by method 1 since the latter produces discrepan- 
cies of 0.49 percent at PFO and 0.33 percent at La Jolla for Mi. 0^ is little affected by tKeans 
other than the north Pacific and has therefore only been calculated using methods 1 and 3. 

A more revealing comparison of the small differences between observation and theory is 
made by comparing only the ocean load part of the tide. Figure 5.3. In the lower part of this 
figure the small load tide part of the total vectors (the boxed region in Figure 5.3) is illustrated 
m greater detail. The origin is therefore the end point of a vector representing the theoretical 
tide of the solid earth with gravimeter factor 1.160 and zero phase lag The arrows represent 
the iK-ean load vectors calculated by method 1 for 0\ and methcxl 2 for Mi. 

In general the agreement between observation and theory is excellent lire calculated 
phases of the O] load vectors are within the scatter of the data with a maximum deviation of 5° 
lor the 54 day sections and T for the complete runs 1. 2. and 3. The calculated amplitudes of 
the load vectors fall outside the scatter of the data points and appear to be approximately 5 
percent too small. For Mi load vectors at both locations the maximum deviation between the 
calculated and observed phase is 8° for the 54 day sections and 3° for the full length runs. At 
PFO Flat both the Mi calculated phase and amplitude are within the data scatter. However, at 
La Jolla the phase is outside the scatter (but only by 3'') and the amplitude is 5 to 10 percent 
too small This result suggests that the ocean tide model may be inadequate for Mi. at least m 
part, because of some variability in its distribution in the neighborhmxl of La Jolla. 
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5.2 Strain Tidest 

Although it is common in earth tide investigations to give tidal responses at different fre- 
quencies, comparison of these with the ocean tides is rarely attempted. Lambert (1974) sug- 
gested that gravity tides measured in Nova Scotia showed the effect of the ocean tide resonance 
in the Bay of Fundy. Warburton and Goodkind (1978), in a study of gravity tides at PFO, 
corrected their measurements for ocean loading by using frequency dependent scaling factors 
for tides in the north Pacific, which were deduced by Beaumont from coastal tide data; they 
found that employing such corrections reduced the apparent scatter in their determinations of 
the tidal response. For such simple scaling to work the load tides must be dominated by the 
effects of a few areas. The Green function for strain falls off with distance more rapidly than 
does that for gravity, so that strain tides recorded near a coast will be mostly affected by nearby 
tides while gravity measurements will not. However, records of earth strain have rarely been of 
sufficient quality or length to enable the tidal response to be determined at any frequencies but 
those of the largest diurnal and semidiurnal tide. The only expectation is the work of Levine 
(1978), who found that the strain tide response varied signifkantly across the diurnal tidal 
band, but only slightly across the semidiurnal band. Some of the diurnal variation could be 
explained by the nearly diurnal core resonance (Wahr, 1979):, the rest presumably showed the 
frequency dependence of the ocean load. No attempt was made to explain this variation in 
terms of ocean tides; because the measurements were made at a mid-continental site, the load 
was not caused primarily by ocean tides in one area. The strain data from PFO were taken near 
the coast, so that ocean load signal is dominated by the load from nearby areas, and a scaling of 
the ocean load vector should have some chance of success. 

In any kind of analysis one always wants to maximize the signal-to-noise ratio (SNR). 
The apparent power spectral density of a sinusoidal signal of amplitude A is 1/2.4^ T, where T 
is the record length. Provided that the noise is stationary, the longest possible record will max- 
imize the SNR. If the noise is not stationary it may be better to eliminate the noisy parts of the 
record, the loss of record length being more than compensated for by the decrease in noise 
level. Visual examination of the strain data showed enough variability in the noise level to 
make a further examination worthwhile. 

In order to reduce the bulk of data to be processed it was decimated from an interval of 6 
to 30 minutes, simultaneously lowpassing it with a filter designed (using the program of 
McClellan et ai, 1973) to have negligible ripple in the tidal bands. To analyze the resulting 
series it was divided into 100 day sections, overlapped so that each one started 800 hours after 
the next. For the estimate of the noise in the frequency band between the diurnal and semidi- 
urnal tides not to be biased by sideband energy from the tides themselves, the data sections 
must be windowed before Fourier transforming (for a discussion of windowing see Harris, 
1976). The 47 t prolate spheroidal window of Thomson (1977) was used, which suppresses side- 
bands to be more than 100 dB below the main lobe. Because this window nearly eliminates 
data from the ends of the series, a 30 percent overlap can be used without the resulting esti- 
mates being highly correlated. A 100 day record gave enough separation between the tides that 
about 60 points in the power spectrum could be used to estimate the intertidal noise level. In 
order to extrapolate this estimate to tidal frequencies the noise power density to be proportional 
to /“‘ was assumed. Though for long periods the intertidal noise level was reasonably con- 
stant, all three strainmeters showed times of high noise. Usually these episodes were at times 
of high drift caused by heavy rain; if the drift rate is high enough it will produce significant 
energy at tidal frequencies. Because the north-south instrument is most sensitive to rainfall it 
showed the greatest number of high-noise episodes. The northwest instrument is the least sen- 
sitive, and only 4 days of data had to be eliminated because of raiti-induced drift, though 47 
days were removed because of high noise produced by malfunctioning electronics. The east- 
west strainmeter showed much higher noise levels for the first year examined than for those 
following, also presumably a reflection of electronics problems. 

+This seciKin is taken from Agnew 0979). 
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The total length of the data series from the north-south and northwest strainmeters was 
2091 4 days, and from the east-west 1926.8 days. All the series began at the same time (day 61 
of 1973); the east-west series is shorter because it was very incomplete for the last 160 days of 
the period analyzed. Because of high noise levels, 46 days were removed from the north-south 
series, 400 from the east -west, and 49 from the northwest. 

Figures 5.4 through 5.6 show power spectra for the three instruments. All three spectra 
display a roughly similar level of background noise falling off with increasing frequency as f \ 
Superimposed on this background are the expected strain tides and also a series of apparent 
tides at exactly 3, 4 . . . cycles per day (cpd). These presumably reflect temperature effects 
manifesting themselves as apparent strains through one of the mechanisms described in Agnew 
(1979). These "thermal tides" are greatest on the north-south instrument and least on the 
northwest, no doubt because the end buildings and part of the vacuum pipe of the latter are 
below ground and thus shielded from daily temperature changes. Close examination of the 
spectra also shows that the level of background noise is about 3 dB less on the northwest 
instrument than on the other two, suggesting that thermal contamination extends outside the 
tidal bands. Unfortunately temperature records are not available for the period of the strain 
data. Levine (1978) found that the strain and air pressure were coherent at periods of a week 
and longer. Cross-spectral analysis of the strain with a 1.5 year record of air pressure showed 
no coherence except at multiples of 1 cpd; this is probably a spurious correlation, owing to both 
the temperature tides and the atmospheric tides being caused by solar heating. 

As part of their study of strain tides recorded in the United States, Beaumont and Berger 
(1975) computed the ocean loading effects at PFO for the 0\ and tides. We have used a 
slightly improved version of their computations provided by Dr C. Beaumont (personal com- 
munication). The ocean tide models used for the Pacific Ocean off the western coast of the 
United States were those derived by Munk et al. (1970) from deep-sea measurements. For the 
rest of the world ocean the A/j cotidal map of Hendershott (1972) and the Oj map of Tiron et 
al. (1967) were used. The Mi tide in the Gulf of California was modeled after the results of 
Hendershott and Speranza (1971). 

Figures 5,7 through 5.9 show how well the computed and observed tidal admittances 
agree. In each figure, the upper part shows the admittance for the diurnal tides, the lower that 
for the semidiurnal tides. The complex admittances are shown as on an Argand plot; positive 
angles (which correspond to phase leads) are plotted counterclockwise. The phases are plotted 
relative to that of the tidal potential, which is the usual reference in ocean tide studies. The 
shear strain is not in phase with the potential; because all three strainmeters partly respond 
to it none of the elastic strains (denoted by E) have zero phase. The computed ocean loads are 
given in actual strains for a Mi (or 0\) tide in the ocean; to normalize these we divided by the 
values of the elastic strain tides, as given by a harmonic analysis of the theoretical series. We 
also corrected the computed load tides for the variation in harmonic amplitudes caused by slow 
changes in the orbital inclination of the Moon. In both figures the vector NP denotes the load 
tide from the area covered by Munk et al. (1970), GC (where present) that from the Gulf of 
California, and RW that from the rest of the ocean. The sum of these is the computed admit- 
tance. shown as an open square. The admittance of the observed tide is shown by a solid dot 
The numbers beside these symbols are the tidal group number, the second element in the 
Doodson number. The 0\ tide is thus marked as -1 in the diurnal plots, and the Mi tide as 0 
in the semidiurnal plots. 

It is clear that the observed and computed tides do not agree; the misfit is more apparent 
than It was in the results of Berger and Beaumont (1976) because the errors on the observed 
admittance are smaller. Berger and Beaumont discussed a number of possible sources of error, 
and concluded that the most likely ones were inaccura ies in the ocean tide models used and 
inadequacies in the representation of strain distortions by inhomogeneities. Though the total 
load from the oceans other than the north Pacific is small, it is made up of a number of larger 
and partly canceling loads, some of which, in particular that from the southeast Pacific, are 
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poorly knc'wn. It is however very I'kely that the corrections for inhomotjeneities are not in fact 
correct because of our ignoiancc of the geology underlying Pinyon Rat. 

The semidiurnal plots for all three strain components, and the diurnal plot for the 
northwest strain, show the admittance for all groups from -2 through 2. The disagreemeii^s 
between observed and computed A/j strain tides suggest that an elaborate modeling effort to 
explain the observed frequency variation would not be warranted. We have instead scaled the 
NP vector by the admittance of the ocean tide at La Jolla (relative to A/j and 0{,, as given by 
Cartwright ei al. (1969), and have scaled the GC vector using some approximate factors derived 
from the results of Stock (1976). Table 5.1 gives the scaling factors used. The resulting com- 
puted strain tides are plotted as open squares, identified by group number. For clarity, the 
observed and computed points have been connected by lines. 

Considering the crudeness of the method, the agreement between the frequency depen- 
dence of the observed and computed admittances is surprisingly good. The results of Munk et 
cl. (1970) show that the frequency dependence of the tides at several points some distance 
f rom the coast is fairly similar to that at La Jolla. As the size of the NP vector shows, the loads 
from this region dominate those of the other areas. The strain tide results thus confirm that 
over this whole region the frequency dependence of the tidal admittance is well represented by 
that at La Jolla, though the growing mismatch between the curves for the tidal groups with fre- 
quencies above A/^ .r jggcsts that this is Icsc true v'ith increasing frequency. 

5.3 Elastic Structure • Local Departures from Radial Symmetry 

Since the strain tides are not in part composed of a direct effect, fractionally they are 
much more sensitive to local clastic conditions than arc the gravity tides. Further because the 
strain Green's function for surface mass loads drops of like away from the source, while the 
gravity Green function attenuates like strain is relatively more sensitive to "nearby" surface 
mass load. 

Before we can deduce much about these surface mass loads, namely, the ocean tides, we 
must correct the observed strain tides for local departures from radial symmetry— a condition 
used to calculate the loading Green’s Functions (Farrell, 1970). 

The distortions of the strain field due to; a) topographic fluctuations of the free surface in 
the vicinity of the strain observatory; b) the cavities in which underground instiut-nems arc 
placed, and c) the lateral variations of elastic parameters caus«l by geological structure, are all 
manifestations of may simply be termed as departures from radiai symmetry. Perhaps the 
most obvious form of lateral heterogeneities are due to geology, but recently the other effects 
have begun to be lully appreciated (Harrison. 1976). The topographic effects are the most 
amenable to model studies simply because the topography is well known and readily available in 
the form of large scale topographic maps. V'ault geometry in the form of three-dimensional 
mine workings is. in general, too complex to model practically Our studies indicate, 
however. that for most underground strainmeter installations the cavity effect is small. Geologi- 
cal effects arc perhaps the most difficult to estimate simply because the three-dimensional distri- 
bution of elastic properties is not well known These effects, unfortunately, arc not necessarily 
small as we will demonstrate of several simple models 

We refer to the “homogeneous applied" strain as the strain applied by the sum of the solid 
earth tide and the (vean load tide to a radially symmetric earth model Conveniently, it can be 
thought of as the strain input to the base of the crust above which lateral variations are 
confined We use the term "inhomogeneous" to describe the strains that result on the free sur- 
face of an irregular heterogeneous earth and these quantities will be denoted by the subscript 
/ 

Finite clement models of lateral variations in the elastic properties of the crust and topog- 
raphy carno accurately represent the total assc.nblage of such variations on the earth This 
docs not. however, prevent realistic modeling because it can be shown by analogy to Sami- 
Venant's principle that each mhomogencity perturbs the strain field only within a surrounding 
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Table 5.1 

Scaling Factors for Ocean Tides 



Gulf of California 

Northeast Pacific 

Group 

Amplitude 

Phase 

Amplitude 

Phase 



(deg) 


(deg) 

1 -2 



1.014 

5.3 

1 -1 



1.0 

0.0 

1 0 



0.850 

-8.5 

1 1 



1.111 

-14.8 

1 2 



1.231 

-20.9 

2-2 

1.0 

-14. 

1.233 

50.8 

2-1 

1.0 

-7. 

1.250 

21.5 

2 0 

1.0 

0. 

1.0 

0.0 

2 1 

1.15 

7. 

0.714 

-17.0 

2 2 

1.30 

14. 

0.908 

11.0 
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region having the same characteristic dimension as the inhomogeneity. At large distances from 
the inhomogeneity the strain field becomes uniform. Furthermore, the length scale of the 
homogeneous tide is, in general, much larger than the length scale of the inhomogeneities. 
Consequently, a valid model is one of a region that surrounds the strainmeter subject to uni- 
form stress, implying a uniform strain field in the absence of inhomogeneities. The region 
should be large enough to encompass all inhomogeneities that will modify the strain field at the 
strainmeter site. This model may not be valid where the strainmeters are sited very close to an 
ocean because local tidal loading produce a homogeneous strain field of wavelength less than or 
equal to that of the inhomogeneities. 

A similar scale argument permits large scale structure to be separated from small scale 
structure in the form of a series of models. For example, the problem of a strainmeter in a 
mine tunnel within a large hill can be separated into two superimposable models; the large scale 
hill, and the smaller scale tunnel. This separation is valid if it can be shown that the stress field 
on each of the sub-models is uniform. This approach is essential when finite element tech- 
niques are used because models that are sufficiently large to include large scale structure yet 
sufficiently detailed to include small scale structure have an unmanageably large number of 
finite elements. 

In the estimation of the topographic effects on the local strain field we consider the strain 
components on an irregular free surface. Thus the strain tensor is of rank one and its com- 
ponents may be represented vectorially. We will refer to the complex strain components in the 
East-North (.v,, axis system (see Figure 5.101 as the vector ?, where 

c, - e,, 

-= P22 

The real part of these complex quantities is the part that is in phase with the local equilibrium 
potential, the imaginary part is that lagging the equilibrium potential by 90°. 

The extensions along the strainmeter directions will be referred to by the vector e" where 
the component €) is the extension in the direction 9\. 

We subscript the quantities ? and ?" with H or I lo distinguish the homogeneous and 
inhomogeneous quantities. The results of the calculations will be ?/, the topographically per- 
turbed applied extensions. 

The results of the finite element modeling of the topography at a particular station are 
expressed by the transformation matrices /' which express the perturbation of the strain com- 
ponents Cif. Thus; 

e, - re^ 

Strain components in the unprimed coordinate system are given by 

e “ 

where R is the rotation matrix 

cos^ <t> 2 s\n<t» cos<t> sin^d» 

—Sind) cosdi cos^<^ - sin‘d> sin d> cos<^ 
sin^d> —2 sin<fi cosd> cos^d> 


Then 




where 

T-- RTR ' 
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Coordinate System 


figure 5 10 Coordinjiie system tor topographic effect uiculaiions 
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These mairices are given in Table 5.2 for PFO. For the cases where <t> is 0® (PFO) the 
rotation matrix is simply the identity matrix. For a station where there were no topographic 
efifects the matrix T wojld also be the identity matrix. 

For PFO we show in Figure 5.11 how various elements of the matrix T actually vary in 
space. This can give some indication when compared with to local topography of good and bad 
sites as far as these effects are concerned. 

A third matrix transformation relates the extensions T to the components of strain 7. 
That is 


where 


C 


r ^ j 

cos'^ <^i 2 sin 9i cos 9] sin'^ 9i 
cos^ 9j 2 sin 9j cos 92 sin^ ^2 
cos^ 92 2 sin 9i cos 0? sin^ 03 


Thus we have the expression 

f, - CRTR '<r '7 „ 

= crc '7„ 

For graphical clarity in the presentation of the final tidal phasor plots (Figure 5.12) we 
express the topographic perturbation effects as an additive vector. 

= ^f/ + T, - 6// 

+ (C7T ‘ - /)?w 

where the additive perturbation vector P is obtained by transformation on the vector?//. 

5.4 Investigating the Temporal Stability 

The third aspect of our earth tide studies directs itself to temporal changes in local elastic 
properties of the crust. 

The observation of changes in the ratio of seismic velocities Fp/ before earthquakes led 
to the suggestion that in the epicentral region large changes in bulk elastic properties may take 
place. A theory is that F^/ F, anomaly, which appears to be predominantly a change in F^. is 
caused by dilatancy in the earthquake source region (Scholz et ai, 1973). Dilatancy, as 
observed in laboratory experiments is an inelastic volume increase of the rock under a stress 
near its failure point. This volumetric strain is due to the formation of micro-cracks. 

We have numerically investigated the effects on the earth’s tidal response due to the pres- 
ence of a zone where F„/ F, changes. We conclude that the earth’s response to tilt and strain 
will be significantly altered by the presence of this anomalous region. Further, because the 
fractional change in the strain and tilt tide response is greater than the corresponding change in 
VJ F,. measurements of tilt and strain tides may provide a more sensitive indicator of dilatancy 
than measurements of seismic velocities. 

Recent Russian results have indicated variations in the strain tide prior to an earthquake 
in Southern Tien Shan of magnitude of 5 a! a distance of 20 km. PFO is ideally suited to such 
a study as it lies between the highly active San Jacinto fault and the main trace of the San 
Andreas fault in a region of high seismicity 

We have used an axisymmetric model of the earth to determine the changes in the 
North-South components of the earth tide in the vicinity of an anomalous region. The dimen- 
sions of the modeled zone, 200 km deep by 2200 km (20®) wide along the earth’s surface are 
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Table 5.2 

Topographic Transformation Matrices 


Station 


T,v 


T,y 

9 


7-,, 

1.G35 

+0.000 

0.005 


Gimp Elliott 

T2. 

+0.000 

1.013 

0.003 

-26° 


hr 

-0.003 

0.005 

1.003 



hr 

0.815 

-0.053 

-0.053 


Pifton Flat 

hr 

-0.053 

0.923 

-0.058 

0° 


hr 

0.003 

-0.063 

1.023 



hr 

0.788 


-0.053 


Pinon Flat 

hr 

0.002 


0.020 

-22° 


hr 

0.010 


1.025 
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Figure 5 12 PFO (Pirton) M 2 and Oj phasors. 
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chosen so that its boundaries are sufficiently distant from the anomalous region that its pres- 
ence docs not effect the strains at the boundaries. We use the boundary conditions appropriate 
to the tide on a laterally homogeneous earth. The results are strictly applicable only to the 
M 2 tide at latitudes 20° to 40°, however, at other latitudes and other tidal components the 
results are not expected to be altered significantly. 

We have used the velocity variations to model Young's modulus, £, and Poisson's ratio, 
cr, within the anomalous region according to the relations 




keeping the density, p, constant. 

The dimensions of the anomalous region are chosen to be 20 km by 40 km but to a cer- 
tain degree, the results scale with these dimensions. Further, we have determined that the 
exact shape and the sharpness of the boundaries do not alter the results greatly. 

The finite element solution provides the radial and tangential displaceme.its of each of the 
surface nodal points. These displacements are differentiated numerically by cubic spline inter- 
polation, to give the north-south components of tilt and strain. Thus, the tilt used here is just 
the slope of the deformed surface. The earth tide tilt measured by a tiltmeter is related to the 
gravitational potential by D, the diminishing factor. £ — where h and k are the 

Love numbers. Changes in Vpf will modify only the h part of the diminishing factor, that is, 
the slope of the Earth’s surface. There will be no change in the direct acceleration and virtually 
no change in k as it is dependent on the deformation of the Earth as a whole and not on the 
small perturbations of the deformation in the anomalous region. 

Figure 5.13 illustrates the results for the radial displacement and north-south components 
of the strain and surface tilt as a function of Vp reduction in the anomalous region. The results 
are normalized by the radial displacement, tilt and strain for the corresponding homogeneous 
model, also shown on the figure. The normalized plots represent the fractional changes in the 
tidal amplitudes caused by the Vp/ Vs anomaly. The surface displacement anomaly is too small 
to be detected with a gravimeter, but both the tilt and strain anomalies will be readily detect- 
able, as we shall demonstrate. 

That the tilt and strain anomalies are both edge effects is immediately obvious. As we 
proceed into a dilaiant region both anomalies decrease,and become undetectable toward the 
center of a very large homogeneous dilatant region. The strain results are particularly 
encouraging, suggesting that a 15 percent decrease in Vp produces a 2 percent tidal strain ano- 
maly even at a distance of 120 km from the edge of the dilatant region. The narrower tilt ano- 
maly is greater than 2 percent at distances less than 30 km from the boundary. The effects of 
smoothing the geometry, burial, or grading the velocity changes gradually are similar and result 
in the removal of the high spatial frequency content of the anomaly. Both the tilt and strain 
anomalies remain significant, however, several tens of kilometers from the boundaries. 

The numerical results may be parameterized in order to relate the changes in elastic pro- 
perties, due to arbitrary Vp and Lj changes, to the surface displacement, tilt, and strain 
anomalies calculated by the finite element models. From the results of this parametrization, we 
conclude that the anomaly maxima are approximately functions of one variable, the Poisson's 
ratio contrast, for tilt and elevation, and the (X ■»- /a)"' contrast for tangential displacement and 
strain. Contrast means the ratio of the values of the parameter inside and outside the 
anomalous region, at the same depth. (X 's the areal strain equivalent of the bulk 
modulus (X -t- 2/3/i), and is termed the "areal bulk modulus." For a 40 km x 20 km anomalous 
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SURFACE DISTANCE (Km) 


f ijiurc ' 13 North south wmponents of the radial displacement, tilt and strain tide amplitude for the model shown at 
the hottimi of the figure Model properties below 26km follow the Gulenburg-Bullen A model Normalizing 
lumiions arc those for the same model without a dilatant inclusion The normalized functions represent the frac 
iional change in radial displacement, tilt and strain, as a function of position, when the P-veloaiy in the dilatant 
region is reduced by 5, 10, 13 and 20 percent 
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region, we can approximately express the maximum normalized tilt anomaly as a function of 
the Poisson’s ratio contrast Aa by 

r„„--1.25Ao- + 2.25 


so that when for example Poisson's ratio changes to half its normal value in the anomalous 
region, the maximum normalized tilt change is approximately 1.6. This means that at the edge 
of such a region, the North-South component of the Mi tilt tide svill be increased by 60 per- 
cent. 


For the same model we calculate that the maximum normalized strain amplitude is 
related to the reciprocal of the areal bulk modulus, by the expression 


*m«x 


0.8A 



-t-0.2 


Thus, for example, if the areal bulk modulus changes by 50 percent across the boundary, the 
maximum normalized strain will be 1.8. niis means the north-south component of the strain 
will be 80 percent bigger than normal 

On the basis of these calculations, Agnew (1979) has examined the strain data from the 
three component laser strainmeier ai the PFO in an attempt either to discover significant 
changes in the tidal admittance that would be caused by a nearby anomalous region, or to set 
practical limits on the detectability of changes in the Mi admittance. 

The admittance for the m-th section analyzed was computed as 


where the sum was over harmonics. The section length, (1093 terms), tapering (Hanning), and 
harmonics averaged over (3, centered on Mi) were the same as were used in computing the 
cross-spectrum. Areas of high noise were not skipped over; as an estimate of the noise we 
used the average power spectral density from 1.70 to 1.78 cpd, which is the high-frequency 
edge of the intertidal band.. The standard error of the admittance estimate is then given 
approximately by 

(Py/n'- 

A 


where T is the record length, Py the noise level, and A the amplitude of the harmonic. 

The phase and amplitude of the Mi admittance, for all the strainmeters. is plotted against 
time in Figure 5.14. There docs not seem to have been any significant departure from the 
mean values for the years 1973-1980. Though the limits on the admittance vary from section 
to section, a change of more than 2 percent in the admittance for 6 months would probably be 
detectable 

Warburton and Goodkind (1978) found that the A/i gravity tide admittance at PFO varied 
by about 1 part in 10 ’ from month to month. The ocean load is about 30 percent of the 
northwest strain tide, and only 2 percent of the gravity tide. If the changes in the gravity tide 
reflected variations in the ocean load they might therefore be expected to alter the strain signal 
by perhaps 15 percent. Unfortunately the noise level is too high to rule out variations of this 
size, though any larger variations would conflict with the strain data 
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